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ABSTRACT 

Aims. Approximately 20% of all spiral galaxies display starburst activity in nuclear rings of a few hundred parsecs in diameter. It 
is our main aim to investigate how the starburst ignites and propagates within the ring, leading to the formation of massive stellar 
clusters. 

Methods. We observed the ring galaxy NGC 7552 with the mid-infrared (MIR) instrument VISIR at an angular resolution of 0'.'3 - 074 
and with the near-infrared (NIR) integral-field spectrograph SINFONI on the VLT, and complement these observations with data from 
ISO and Spitzer. 

Results. The starburst ring is clearly detected at MIR wavelengths at the location of the dust-extincted, dark ring seen in HST 
observations. This "ring", however, is a rather complex annular region of more than 100 parsec width. We find a large fraction of 
diffuse [Ne n] and PAH emission in the central region that is not associated with the MIR peaks on spatial scales of ~ 30 pc. We do not 
detect MIR emission from the nucleus of NGC 7552, which is very prominent at optical and NIR continuum wavelengths. However, 
we have identified nine unresolved MIR peaks within the ring. The average extinction of these peaks is Ay = 7.4 and their total 
infrared luminosity is L m = 2.1 x 10 10 L G . The properties of these peaks are typical for MIR-selected massive clusters found in other 
galaxies. The ages of the MIR-selected clusters are in the range of 5.9 ± 0.3 Myr. The age spread among the clusters of 0.8 Myr is small 
compared to the travel time of ~ 5.6 Myr for half an orbit within the starburst ring. We find no strong evidence for a scenario where 
the continuous inflow of gas leads to the ongoing formation of massive clusters at the contact points between galactic bar and starburst 
ring. Instead, it appears more likely that the gas density build up more gradually over larger ring segments, and that the local physical 
conditions govern cluster formation. We note that the fundamental limitation on the accurate derivation of cluster age, mass and IMF 
slope is the lack of higher angular resolution. Resolving the highly embedded, massive clusters requires milli-arcsecond resolution at 
infrared wavelengths, which will be provided by the next generation of instruments on extremely large telescopes (ELTs). 

Key words. ISM: kinematics and dynamics, HII regions - Galaxies: nuclei, starburst, star clusters - Infrared: ISM 



1. Introduction ters, accumulate and concentrate in these ring-like structures. 

Unfortunately, the situation is not clear, partially due to the lack 
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(N ; Massive star formation in galaxies occurs in various appearances of high resom tion observations at longer wavelengths and par 

^ , and modes, from localised H II regions and blue compact dwarf tially due to the comp i ex ity of the dynamics of gas and dust, and 

^ ■ galaxies to luminous starbursts in ultra-luminous infrared galax- hence several scen arios have been discussed in the literature. 

k> '■ ies ( ULIRGs ) and sub-millimeter galaxies. One particularly in- According to IShlosman et all (fl990l) . the torque from the 

^ . teresting starburst mode shows a ring-like morphology around i arge . S cale bar weakens near its centre where the influence of 

H ■ the galactic nucleus (e.g., | Sarzi et al.| | 2007t |Bokeretal. | |2008t me bulge begins t0 dominate the disk pote ntial. The gas concen- 

- - | MazzucaetalJ | 2008|) . These systems are not uncommon - in ^ in ring . like resonance s where th e gas is free of gravita- 

approximately 20% of all spiral galaxies star-formation occurs tional torques In contrast> KenneviLiD dS) explain the rings 

primarily in these rings of typically a few hundred parsecs in afj a conseque nce of a previous nuclear starburst which has used 

diameter ( jKnapen| | 2005|) . These large structures of intense star up the gas in the centre, leaving the nuclear ring as a remnant. 



formation are primarily found in barred spiral galaxies of types Recentlyi |y a isanen et alJ (T20ll reported on an outward propa 



Sa-Sbc. In fact, bars are expected to account for about 3.5 times gating ring of star formation in NGC 1614 resulting from a re l a - 
more triggered cen tral star formation than galaxy interactions t j ve jy ma j or merger event 

^ T ^° ^j l2^4^gPg ma g net j c fields dO^G in the ring Arguably the most common explanation is that the gas ac- 



|| i ,., , rtiguauiy me mosi common exmaiiauon is mai me gas ac- 

ofNGC7552, Beck 201 1) are likely to create magnetic stress , . , ■ • a \ . , Jr . , 

. ' 1 ?' v , , J fe cumulates during the inflow around the radii at which the stel- 



that causes inflow of gas toward the center. lar ^ experie nce dynamical reson ances with the rotating 

The key question is under what conditions the large amounts bar poten tial (e.g., Binn ev & Trem1dn^l2008l) . In the central re- 

of gas and dust, necessary to form the massive young star clus- gion this typica u y happens at the so-called inner Lindblad res- 

onance (ILR), arisin g from the interplay between the bar and 

Send offprint requests to: B. R. Brandl, brandl@strw.leidenuniv.nl the stellar orbits (e.g jTelesco et al.ll993blButa & Combesll996t 
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Bok er et al.l 120 08). The gas is trapped between an inner ILR 
and an outer ILR. If the potential allows only one ILR, the 
gas will concentrate directly at the ILR, if the system has two 
ILRs, the gas will acc umulate at the inner ILR - (see review by 
iButa & Comb"esl[T996l for more details). 

However, more recently iRegan & Teubenl (120031) claimed 
that the common assumption that nuclear rings are related to 
an ILR is incorrect, and showed in simulations that there is no 
resonance at the inner Lindblad resonance in barred galaxies. 
They argue that confusion has developed between the orbit fam- 
ily transition that occurs at these radii in weak b ars and a true 
resonance. Furthermore, Regan & Teuben (2003) found that the 
radius of nuclear rings decreases with time, either because the 
rings accumulate lower specific angular momentum gas or be- 
cause of dissipation at the contact point of the bar dust lane and 
the nuclear ring. 

Obviously, observations may be able to provide further in- 
sights. A very good example of a starburst ring galaxy is 
NGC7552, a nearl y face-on, barred sp iral galaxy at an incli- 
nation of ~ 28° (Feinstein et al.l Il990t) . Originally classified 
as amorphous ( Sersic & Pastoriza ll965l) it appears in the RC3 
catalogue (de Vaucouleurs et al.l 119911) as morphological type 
SB(s)ab and in the RSA cataloeue lSandage & Tammannf l981) 
as type SBbc(s). Both designations indicate a barred spiral (SB) 
galaxy with spiral arms that spring from the ends of a bar 
(s), and an outer ring (R'). Fig. Q] shows a 3-color image of 
NGC 7552 to illustrate the global morphology and the relevant 
scales. Throughout th is paper we assume a distance of 19.5 Mpc 
dTullv & F isher 1988); at this distance, 1" corresponds to 95 pc. 

NGC 7552 has been previousl y classified as a L INER 
galaxy, based o n the [011/16300 line dPurret & Bergeronlll988l) . 
Neither X-rav (|Liu & Breemanl 12005 ) nor NIR observations 
dForbes et al.lll994ah did reveal significant activity from a galac- 
tic nucleus. In f act, models of gas flow in barred galaxies 
dPiner et al . 1995) have shown that the gas flow, necessary to fuel 
an active galactic nucleus, can be interrupted by a nuclear ring. 
The global spectral energy di stribution (SEP), from UV to ra- 
dio, is well characterised (e. gjDale et al. 120071). From the IRAS 
colours of the entire galaxv. lSanders et all (12003) estimated a to- 
tal infrared luminosity of 8.5 x 10 10 L o (adjusted for the assumed 

distance). 

1990t 



possible to study galaxies in the N band (8-13 /mi) at an angular 
resolution of 073 under excellent seeing conditions^. 



Several previous studies (e.g., Feinstein et al 



Forbe s et al.lll994bt ISchinnerer et all 19971: ISiebenmorgen et alj 



2004: Ivan der Werf et al.l 120061) provided strong evidence for a 
starburst ring in the centre of NGC 7552 (section 13.1) . While 
the macroscopic picture of how these rings form seems well 
established, the morphology of the central region and the 
microscopic picture of how, when and where the starburst 
ignites and propagates are rather uncertain and require more and 
better observational constraints. It is the main aim of this paper 
to investigate these details. 

MIR observations pr ovide excellent diagnostics to detect 
embedded cluster s^] (e.g. iGoruan et all 1200 U iBeck et al 2002; 
Vacca et all 120021) and q u antify massive star formation (e.g. 
Genzel & Cesarskvl 120001: iPeeters et all |2004|) . Previous stud- 



ies of NGC 7552 at MIR wavelengths have been severly lim- 
ited in spatial resolution, with angular resolutions o f 0.8" - 1" 
(ISchinnerer et all 119971: ISiebenmorgen et al.l 120041) and ~ 4" 
(Pal e et al .1120071) . The availability of the VLT Spectrometer and 
Imager for the Mid-Infrared (VISIR) at the ESO VLT, makes it 




1 Throughout this paper we use the physical term 'cluster' and the 
observable 'MIR peak' as synonyms. Whether or not this is correct will 
be discussed in section 1431 



Fig. 1. Colour composite o f NGC 7552 made fro m images pro- 
vided by the SINGS team dKennicutt et alj|2003l) in H-a (blue), 
and the B and I bands (green and red, respectively), taken on 
the CTIO 1 .5 meter telescope. North is up and East to the left. A 
scale bar indicating angular and absolute distances is also shown. 



In this paper we complement our mid-IR data with NIR im- 
ages at similar angular resolution. The latter data were taken 
with SINFONI, a NIR integral-field spectrograph at the ESO 
VLT. Combining near- and mid-IR photometry will provide es- 
sential information on both, the stellar sources that provide most 
of the energy output, and their dusty, re-radiating environment. 
The observations presented here, with their unprecedented high 
angular resolution, will allow us to study the morphology of the 
starburst ring, the distribution and properties of its star clusters, 
and their formation history. 

This paper is structured as follows. Section |2] describes the 
observations and data reduction. Section [3] presents the photo- 
metric and spectroscopic results, and sect. |4]provides a detailed 
discussion of our results. Finally, sect. summarises our main 
conclusions. 



2. Observations and Data Reduction 

2.1. VISIR 
2.1.1. Imaging 

The imaging data were acquired with th e VLT Imager and 
Spectrograph for the mid-IR (VISIR, lLagage et all 12004 
iPantin et al]|2005l) . mounted at the Cassegrain focus of the VLT 
Unit Telescope 3, Melipal. The data were obtained in period of 
March 2005 to January 2006 in nine different filters. TableQ]lists 
the relevant filters, their central wavelengths and spectral widths. 
A cross-comparison with the spectro-photometry (section |2.1.2I > 
shows overall very good agreement in eight narrow and broad- 
band filters, but with the exception of the SIV_1 filter, in which 
the measured flux density appears to be about 60% too high (see 
Fig.HJi. However, the 9.73 -9.91/mi SIV_1 filter has intrinsically 
very low transmission, resulting in a four times lower nominal 



2 A seeing of OC'3 at lOyum corresponds to (X'55 optical seeing. 0?3 is 
also the diffraction-limited resolution of an 8 m telescope at 10 //m. 
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Table 1. The spectral properties and the measured flux densities in units of mJy through the nine used VISIR filters. In addition 
to the individual MIR peaks Ml through M9 we list their summed contribution (ZMj), the integrated fluxes within an annulus of 
1'.'5 < r < 4", and the image background noise (<x in units of mJy arcsec -2 ). Upper limits are defined as three times the image noise 
within a circle of 077 in diameter. 



Filter 


SIC 


PAH1 


Arlll 


SIV_1 


SIV 


SIV.2 


NelLl 


Nell 


NeIL2 


A c [fim] 


11.85 


8.59 


8.99 


9.82 


10.49 


10.77 


12.27 


12.80 


13.03 


AA [/mi] 


2.34 


0.49 


0.14 


0.18 


0.16 


0.19 


0.18 


0.21 


0.22 


Ml 


116 ± 13 


67 ± 10 


80+ 11 


119 ± 14 


72 ± 13 


119 ± 14 


184 ± 14 


282 ± 18 


200 ± 15 


M2 


90 ± 12 


51+9 


41 ±9 


54 ± 12 


22 ± 10 


66 ± 12 


111 ± 11 


198 ± 15 


132 ± 12 


M3 


41 ± 10 


<50 


< 59 


<78 


<78 


<76 


49 ±9 


111 ± 12 


62 ±9 


M4 


44 ± 10 


<50 


< 59 


<78 


<78 


<76 


47 ± 8 


115 ± 13 


53 ±9 


M5 


<61 


<50 


< 59 


<78 


<78 


<76 


<44 


54 ± 10 


31 ± 7 


M6 


<61 


<50 


< 59 


<78 


<78 


<76 


40 ±8 


63 ± 10 


34 ±8 


M7 


70 ± 11 


83 ± 11 


42 ± 9 


<78 


<78 


<76 


90 ± 11 


141 ± 14 


110 ± 11 


M8 


< 61 


<50 


< 59 


<78 


<78 


<76 


28 ±7 


60 ± 10 


38 ±8 


M9 


41 ± 10 


<50 


< 59 


<78 


<78 


<76 


41 ± 8 


78 ± 11 


48 ± 8 


EMi 


402 ± 66 


201 ± 30 


163 ± 29 


173 ± 26 


94 ±23 


185 ± 26 


590 ± 76 


1102 ± 113 


708 ± 87 


Ring 


2023 ± 81 


1336 ±65 


583 ±59 


1114 + 79 


to 


1262 ± 80 


2326 ± 76 


3700 ± 94 


2435 ± 76 


(T 


37 


30 


27 


47 


47 


46 


27 


35 


27 



It was not possible to perform aperture photometry because of insufficient signal-to-noise. 



sensitivity with respect to comparable VISIR filters^, and was 
thus not included in the analysis. 

The data were acquired in service mode queued observa- 
tions under stable weather conditions with an optical seeing of 
< 078 and an air mass < 1.5. We selected the pixel scale of 
07127 pixel 1 , which provides a total FOV of 3273 x3273. The 
observations were performed using the standard chopping and 
nodding technique, which removes the time-variable sky back- 
ground, telescope thermal emission and most of the so-called 1/f 
noise. We used a chopper throw of 30". All observations were 
bracketed by photometric standard star observations. We note 
that the chopping/nodding technique, while effectively reducing 
the background, may also eliminate any large-scale diffuse neb- 
ular emission from the region. VISIR imaging, however, recov- 
ers nearly 50% of the total [Ne n] line flux measured by Spitzer 
(TableO and hence we do not consider chopping/nodding to in- 
troduce a significant source of systematic photometric errors. 

The final images are the result of shifting and combining 
the individual chopping cycles and were obtained using the 
VISIR pipeline (version 1.3.7). Practically all of the final im- 
ages showed some stripes due to high-gain pixels. We removed 
these stripes using the reduction routine destripe (Snijders, pri- 
vate communication), which is written in IDL. The conversion 
from counts to physical flux units was derived from the photo- 
metric standards HD 216032 and HD 218670. The angular reso- 
lution varies between 073 and 074, very close to the diffraction 
limit of VLTat 10 /mi. 

Fig. |2] illustrates the image quality of the resulting maps in 
the continuum subtracted [Ne n][12.8^m] line, the 8.6yum fea- 
ture of polycyclic aromatic hydrocarbons (PAHs) and the MIR 
continuum. Although the PAH map suffers from fixed pattern 
noise, the [Ne n] and continuum maps illustrate the significant 
impro vement in comparison to previous mid-IR observations 
(e.g., ISchinnerer etal.ll 19971 Fig. 4). 

The starburst ring and several individual knots are clearly 
visible in the 12.8^m filter image in Fig. [3]a. We have identi- 
fied nine unresolved peaks, which we label Ml to M9. The posi- 



3 http://www.eso.org/sci/facilities/paranal/instruments/visir/inst/ 




Fig. 2. Left: Continuum subtracted [Nen][12.8/mi] line image. 
Center: 8.6/im PAH image Right: 12yum continuum images 
taken in the SiC filter. The PAH image suffers from fixed pat- 
tern noise, but shows one bright peak, M7. All three images are 
displayed in linear intensity scaling; North is up and East is to 
the left. The bar indicates a scale of 2" or 190 pc. 

Table 2. Absolute and relative coordinates of th e MIR peaks 
with re spect to the radio nucleus (Nuc) from [Forbes et al.l 
(fl994bh . 



Source 


a (J2000) 


6 (J2000) 


Aa(") 


ash 


Nuc 


23*16"' 10571 


-42deg35'04."4 


+0.0 


+0.0 


Ml 


23'' 16"' 10574 


-42deg35'01."6 


+0.3 


+2.8 


M2 


23'' 16 m 10567 


-42deg35'01."8 


-0.4 


+2.6 


M3 


23'' 16"' 10557 


-42deg35'02."6 


-1.6 


+ 1.8 


M4 


23* 16 m 10/56 


-42deg35'01."8 


-1.7 


+2.6 


M5 


23 A 16"' 10548 


-42deg35'02."3 


-2.6 


+2.1 


M6 


23* 16"' 10549 


-42deg35'06."l 


-2.5 


-1.7 


M7 


23'' 16 m 10585 


-42deg35'06."7 


+ 1.6 


-2.3 


M8 


23'' 16"' 10594 


-42deg35'04."8 


+2.5 


-0.4 


M9 


23'' 16 m 10584 


-42deg35'03."4 


+ 1.5 


+ 1.0 



tions of these peaks are listed in Table [2] For each VISIR filter, 
the fluxes of each knot were measured performing aperture pho- 
tometry within circular apertures with a diameter of 077, about 
twice the resolution limit. The photometric fluxes are all tabu- 
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lated in Table Q] Ml and M2 are the brightest knots, clearly de- 
tected in all the VISIR filters. The quoted uncertainties are those 
derived from the aperture photometry. The summed contribu- 
tion of all the individual knots (EM;) and the total photometry 
obtained within an 3"-8" annulus that covers the starburst ring 
(section [3~TI ) are also listed. 

The [Nen] 12.8 yi/m line emission of the individual clus- 
ters has been estimated from the 'NeII_l', 'Nell' and 'NellJZ' 
narrow-band filters. We interpolated the continuum level at 
the central wavelength of the 'Nell' filter (12.8yum) from the 
'NelLl' and 'NeIL2' filters. The resulting [Nen] line fluxes are 
listed in Table [3] 

Similarly, we have obtained approximate measurements of 
the 8.6 yum emission feature of PAHs. The three most luminous 
peaks in the 'PAH_1' filter are the mid-IR sources Ml, M2 and 
M7. Here we have only used the 'ArlH' filter to measure the con- 
tinuum level and assumed a flat continuum in F v . The results are 
listed in Table [5] With only one data point to constrain the base- 
line image the fluxes are more uncertain and should be compared 
to spectroscopic measurements, where available. 



J Band Integrated Spectrurr 






Br, 

L .. : 




i v - 




V 


y 



2.14 2.16 
Rest Wavelength [/xm] 



Fig. 5. Segments of the spatially integrated spectra at J (top) and 
K band (bottom) to illustrate the signal-to-noise. The strongest 
emission lines ([Fen], Pa-/J, H2 and Br-y are labelled. 



2.1.2. Spectroscopy 

A low-resolution spectrum at N-band was obtained of the bright- 
est peak Ml on June 16th, 2005. The slit with a width of 0775 
was centred on Ml and oriented along PA. =-10°. The spectral 
resolution is R » 185 - 390. Four spectral settings (/l ce ntre = 8.8, 
9.8, 11.4, and 12.4/vm), overlapping by at least by 15%, were 
used to cover the full N-band. The total on-source integration 
time was 30 minutes for each spectroscopic setting. Since the 
images show no large-scale diffuse emission, chopping and nod- 
ding on-slit were applied with a 12" chopper throw. The early- 
type AO stars HIP 109268 and HIP 1 13963 were observed before 
and after each target observation for spectro-photometric flux 
calibration. 

We reduced the data with the VISIR pipeline, which includes 
subtracting the chopped/nodded pairs, correcting for optical dis- 
tortion, wavelength calibration and extracting the spectra from 
a 1727 aperture in the spatial direction. Absolute flux calibra- 
tion was obtained by integrating the standard star spectra over 
the VISIR filter bands and normalising their fluxes to the ones of 
the photometric standards. The accuracy of the flux calibration is 
approximately 20%. The resulting spectrum is shown in Fig.|4] 

2.2. SINFONI 
2.2.1. Data Analysis 

Observations of NGC 7552 were also taken with the 
Spectrograph for INtegral Field Observations in the Near 
Infrared (SINFONI) at the VLT. SINFONI provides spatial and 
spectral data in the form of a data cube covering the J, H, and 
K bands. The SINFONI instrument is mounted at t he Cassegrain 
focus of the VLT Unit Te lescope 4, Yepun (e.g. iBonnet et alJ 
12004 iGillessen et alj |2005). The observations were made in the 
J, H, and K bands using a spatial resolution setting of 0725 
pixel -1 , corresponding to a field of view of 8"x8", and were ob- 
tained on August 27th, 2005 without adaptive optics correction. 
The average air mass at the time of observation was 1.23 and the 
DIMM seeing between 1" and 173. Each galaxy had an integra- 
tion time of 300 seconds per band and a spectral resolution R of 
2000, 3000, and 4000, for the J, H, and K bands respectively. 

The r eduction of the data was performed with the SINFONI 
pipeline (Modiglia ni et alJl2007h . version 1.7.1. This included 



the removal of detector signatures (geometric distortions, bad 
pixels, pixel gain variations, etc.), the sky emission correction, 
the wavelength calibration and the image reconstruction from 
the image slices. Flux calibration and removal of the contami- 
nation from telluric lines was performed using observations of 
the standard star HIP 025007. In Figure [3]c,d we show the con- 
tinuum subtracted Pa-/3 and Br-y line maps and representative 
sections of the J and K spectra to indicate the signal-to-noise. 
Although the H band was included in our observations, it is suf- 
fering heavily from OH atmospheric line contamination and not 
relevant for the following analysis. 

We have also constucted the K band continuum image from 
integrating the SINFONI spectra across the K band without in- 
cluding the emission from strong lines. The result is shown in 
Fig- 0b and discussed in section [37X1 

Table [3] lists the derived Pa/3 and Br-y line fluxes, and the 
equivalent widths (EWs) of Br-y, averaged over circular aper- 
tures of 077 centred on each mid-IR peak. The integrated line 
fluxes within the 3"-8" annulus are also given in Table [3] 
We emphasize that our measured Br-y line flux of 6.55 x 
10 -21 WcirT 2 , integrated over the central region, is in excellent 
agreement with the corresponding values fr om previous studies 
by di fferent authors: 6.5 x 10™ 21 WcirT 2 dMoorwood & Olival 
1990), 6.4 x 10~ 21 WcirT 2 dForbes et all Il994ah . and 5.7 x 



10" 21 Wcnr 2 dSchinnerer et al.lll997l) . 



2.2.2. Uncertainties 

We estimate the uncertainty in the measured Pa-/? and Br-y line 
fluxes to be 15% (Table [3j. This uncertainty is a combination 
of the typically 10% absolute uncertainty on the SINFONI line 
fluxes - based on the residual slope in the standard star spectra 
after the division by the appropriate black body - and the sys- 
tematic uncertainty in the peak position (estimated from shifting 
the extraction apertures by +1 spaxel). 

The uncertainty of the EW of Br-y is only 5% for two rea- 
sons. First, because an EW is a ratio of two measurements at 
similar wavelenths, the absolute photometric errors cancel out. 
Second, because the EW changes (with respect to the line flux) 
less rapidly with distance from the position of the IR peaks the 
EW is less sensitive to the positioning of the extraction aperture. 
The EW of Br-y is mainly of importance for the derivation of 
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2 

id 
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U 

cu 
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u 

< 

-1 
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4 3 2 1 -1-2 -3 4 3 2 10 -1-2 -3 -4 




4 3 2 10 -1-2 -3 -4 4 3 210 -1-2 -3 -4 

Arc Seconds Arc Seconds 



Fig. 3. a) upper left: VISIR image taken in the 12.8yimi filter. The 1 - <x RMS is approximately 0.57 mJy pixel 1 . The plus signs 
indicate the locations of the MIR peaks Ml through M9, discussed in this paper, b) upper right: K-band continuum image of 
NGC 7552, composed from the reconstructed SINFONI data cube, c) lower left: Reconstructed and continuum subtracted SINFONI 
Pa-/? image, d) lower right: Reconstructed and continuum subtracted SINFONI Br-y image. All four figures indicate the image scale 
in arcseconds (1" corresponds to 95 pc) with respect to the radio continuum center (see Tableland section |2~4] |. 



the cluster ages . We note that all of the measurement uncertain- 
ties are small in comparison to the systematic errors given by the 
model assumptions (see sections [331 and |4~3T > . 

2.3. Spitzer-IRAC 

We have also analysed archival IRAC/Spitzer images of 
NGC 7552 from the SINGS Legacy Science Program (e.g. 
iDaleet all 120051) . The images are shown in Fig. [6] Due to the 



smaller telescope aperture the angular resolution is much lower, 
approximately ~2'.'4. Nevertheless, two distinct, bright peaks in 
the ring are clearly visible. We have also performed aperture 
photometry using the standard IRAC 10 pixel (12'.'2) aperture 
radius with a sky annulus of 10 - 20 pixels. According to the 
IRAC Data Handbook (version 3.0) the photometric errors are 
dominated by the absolute calibration uncertainty of 10% result- 
ing from the IRAC filter bandpass resp onses in sub-array mode 
dOuijada et alj|2004c iReach et alj |2005). The derived photomet- 
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2 4 6 8 10 12 14 16 8 9 10 11 12 13 

Wavelength (yum) Wavelength {/J.m) 

Fig. 4. Left: The mid-IR SED of the total integrated flux within the 1 ."5 < r < 4" ring (black solid line) and the sum of our nine 
peaks EM; (black dashed line). The integrated flux of the entire galaxy as measured through the SiC filter (Table [1} is indicated by 
the d iamond. The filled circ les show the IRAC photometry using a 12'.'2 aperture. The filled square is the 12 /mi IRAS flux taken 
from Gil de Paz et"al] J20 06). The spectra from ISOPHOT (blue line) and Spitzer-IRS (red line) were taken through apertures of size 
24" x 24" and 27" x 43", respectively. Right: The VISIR low-resolution N band spectrum of peak Ml (black line) in comparison 
to the Spitzer-IRS spectrum (red line). The latter covers the central 27" x 43" of NGC 7552 and was normalised to the former at 
10//m (scaled down by a factor of 13). In both spectra, the [Ne 11] fine-structure line at 12.8 pm is clearly present but the strength of 
the PAH 1 1.3 /mi band differs significantly. In addition, the blue crosses indicate the VISIR narrowband photometry of Ml and the 
black diamond corresponds to the broadband SiC filter measurement. 

Table 3. The physical properties of the clusters Ml - M9 derived from O'.'l apertures centered on the mid-IR peaks: Pa/3, Bry (both 
measured and extinction corrected), equivalent widths of Bry, and the [Nen] line, and 12 //m continuum (NelLl filter) fluxes. The 
derived values of extinction Ay, number of Lyman-continuum photons A^ vc , the equivalent number of 07V stars, the total infrared 
luminosity, and the approximate cluster ages are also listed. For a discussion of the quoted uncertainties see section [2] 



Source 


PayS 


Bry 


(BryW 


EW(Bry) 


[Nen] 


12 yum cont. 


A v 


N Lyc 


Noiv 




Age 






(10- 23 Wcnr 2 ) 




(A) 


(10- 21 Wcm 


2 ) (mJy) 


(mag) 


(10 51 




(xlO 3 ) 


(10 9 L o ) 


(Myr) 


Ml 


27.1 ±4.1 


12.4 ± 1.9 


23.5 ±4.9 


65.5 ± 3.3 


25 ±9 


184 ± 14 


8.1 


7.8 ± 


1.6 


1.4 ±0.3 


6.0 


5.6 


M2 


22.1 ±3.3 


9.6 ± 1.4 


17.6 ±3.7 


45.8 ±2.3 


21 ±7 


111 ± 11 


7.6 


5.8 ± 


1.2 


1.0 ±0.2 


4.0 


5.8 


M3 


15.1 ±2.3 


7.7 ± 1.2 


15.4 ±3.2 


42.0 ±2.1 


15 ±5 


49 ±9 


9.0 


5.1 ± 


1.1 


0.9 ± 0.2 


1.9 


5.9 


M4 


22.4 ± 3.4 


9.8 ± 1.5 


17.9 ±3.8 


74.5 ± 3.7 


18 ±5 


47 ±8 


7.7 


5.9 ± 


1.2 


1.1 ±0.2 


1.6 


5.5 


M5 


10.3 ± 1.5 


4.4 ± 0.7 


7.9 ± 1.7 


43.9 ±2.2 


7±4 


< 44 


7.5 


2.6 ± 


0.5 


0.5 ±0.1 


0.9 


5.8 


M6 


11.3 ± 1.7 


4.6 ± 0.7 


8.2 ± 1.7 


15.4 ±0.8 


8±5 


40 ±8 


7.1 


2.7 ± 


0.6 


0.5 ±0.1 


1.0 


6.3 


M7 


17.7 ±2.7 


7.9 ± 1.2 


14.7 ±3.1 


20.4 ± 1.0 


11 ±6 


90+11 


7.9 


4.9 ± 


1.0 


0.9 ± 0.2 


3.3 


6.2 


M8 


13.3 ±2.0 


4.5 ± 0.7 


7.0 ± 1.5 


15.5 ±0.8 


7±5 


28 ±7 


5.6 


2.3 ± 


0.5 


0.4 ±0.1 


1.1 


6.3 


M9 


12.1 ± 1.8 


3.9 ± 0.6 


5.8 ± 1.2 


16.8 ±0.8 


9±5 


41 ±8 


5.2 


1.9 ± 


0.4 


0.3 ±0.1 


1.4 


6.3 


EM; 


151.4 ±22.7 


64.8 ± 9.7 


118.0 ±24.8 




121 ± 17 


590 ± 76 


7.4 


39.1 h 


:8.2 


7.0 ± 1.5 


21.3 


5.9 ±0.3 


Ring 


2020.8 ±303.1 


654.8 ± 98.2 


1002.0 ± 210.4 




382 ± 20 


2326 ± 76 


5.3 













ric fluxes are 198 mJy, 1 85 mJy, 624 mJy and 1 809 mJy for IRAC 
bands 1 through 4, respectively. These flux densities are also in- 
dicated in Fig. [4] 

2.4. Relative Astrometry 

Our analysis is based on multi-wavelength data, and the co- 
alignment and relative astrometry of the various maps is of ut- 
most importance , as illustrated in Fig . [7] Our starting point is 
the radio map of iForbes et al.l (1 199 4b) with source B (radio nu- 
cleus, Tabled as reference point. The absolute astrometry of the 
VISIR images was then derived from shifting the [Ne 11] image 
to best match the radio map. Next, we compared the morphology 
of the Br-y emission (Fig.[3]d) with the [Ne 11] map (Fig.|3]a) and 
found an excellent agreement concerning the cluster peak po- 



sitions. We then shifted the SINFONI maps by a small amount 
(within the VLT pointing uncertainties) to match the Br-y and 
[Nen] maps. 

The adjusted K band continuum image (Fig. 0b) shows the 
K band nucleus at RA = 23 /l 16"W69; 5 = -42deg35'04."8, 
which is about 0."5 to the South from the radio nucleus. We em- 
phasize that the offset between apparent K-band nucleus and ra- 
dio center (see also Fig. [3} is bootstrapped from the Br-y image, 
which is part of the IFU data from which the K-band continuum 
was reconstructed. 

The HST images shown in Fig.|7]were combined on the basis 
of the HST reconstructed pointing, and the [Ne n] 12.8 /mi con- 
tours have been adjusted to the radio map (as described above) 
and then overplotted onto the HST coordinate frame. Like the 
K-band continuum, the HST images show the nucleus slightly 
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offset to the South with respect to the active starburst ring. 
However, even if this offset would indicate the uncertainty in 
the absolute astrometric registration, the scientifically more rele- 
vant uncertainty here is the relative astrometry between emission 
peaks, which is approximately two VISIR pixels or (X'25. 



3. Analysis 

3.1. Appearance of the Ring at Different Wavelengths 

As mentioned in the Introduction, the morphology of the inner 
region of NGC 7552 is quite complex. Based on imaging at B VI 
bands, Feinste in et al. I (fl990l) found a nuclear region of low ex- 
tinction within 4'.'9 in diameter, surrounded by a massive gas 
and dust ring, whic h is responsible for most of the far infrared 
flux. Subsequently, [Forbes et al.l d!994bl) analysed 3 and 6 cm 
radio maps, combined with H-a images and found a starburst 
ring of diameter 9" x 7" (0.86 kpc x 0.67 kpc). Fo r comparison, 
the sample of 22 nuclear rings of Maz zuca et al.l (I2008I) shows 
diameters ranging from 0.4 kpc x 0.4 kpc to 3.4 kpc x 2.0 kpc. 
The ring in NGC 7552 looks alm ost circular in projection, due 
to its inclination of only ~ 28° dFeinstein et al.lll990h . and has 
a wi dth of 3" - A". It c ontains about 60% of the 6 cm radio 
flux. iForbes et al.l (Il994bl) state that this ring cannot be seen in 
single colour optical images largely because of dust extinction. 
High resolution NIR im ages taken with the SHARP 1 camera 
( Schi nnerer et al.llT997h showed an inner ring-like structure with 
a diameter of 7" - 8'.'5, and similar structu r es have been seen at 
MIR wavelengths i Schinnerer et al.l 1 1997t ISiebenmorgen et all 
l2004Uvan der Werf et al.ll2006h . 

In order to assess the correspondence between the morphol- 
ogy seen in our MIR observations and other wavelengths at high- 
est angular resolution, we have created a false colour image from 
archival HST/WFPC2 data in the F336W and F439W broadband 
filters and the H-a (F658N) narrowband filter. Fig.Qshows the 
central region, together with the contours of the VISIR [Ne n] 
12.8 /urn emission. We find an excellent correspondence between 
the heavily dust extincted "dark ring" and the [ Nen] emission. 
The H ST images also confirm the early results of Fei nstein et alj 
( 1990) of a central (< 4'.'9) region of low extinction surrounded 
by a massive ring of gas and dust. However we emphasize two 
findings: 

First, the optically dark ring is the origin of most of the MIR 
emission (Fig. [3]), and hence, the most recent population of mas- 
sive stars in the center of NGC 7552 is still hidden at optical 
wavelengths. We also note that this happens on large scales of 
hundreds of parsecs, not just locally for an individual cluster. 




Fig. 7. Colour composite of the nuclear region of NGC 7552 
made from archival HST/WFPC2 images in the F658N (H-a, 
blue) F336W (green), and F439W (red) filter. The yellow con- 
tours outline the [Ne n] 12.8 /mi emission as detected with VISIR 
(Figf3]a), the scale bar indicates 2" or 190pc. North is up and 
East to the left. 



Furthermore, the Br-y map (Fig.[3]d) resembles the MIR [Nen] 
map (Fig.[3]a) very well. 

Second, the reconstructed K band continuu m map (Fig .[3]b), 
which is very similar to the broadband image of Schi nnerer et al.l 
(I1997L Fig. 2), does not resemble the Br-y map at all, but looks 
smilar to the broadband optical colours seen by HST (Fig. |7). 
The northern K band continuum peaks (which must not be con- 
fused with the MIR peaks Ml & M3, which are at larger radii 
from the nucleus) coincide with the bright regions in the HST 
image. Most of the K band continuum emission comes from 
within the inner 3". 

As far as the structure within the ring is concerned we also 
observe significant changes with wavelength. The IRAC images 
(Fig. [6]) show two peaks, the northern one, which is close to 
our peaks Ml and M2, and the southern one, which coincides 
with the position of M7. The two peaks have similar brightness 
at 8 //m, while the southern one becomes the dominant peak at 
shorter wavelengths. Our strongest 12 ytzm source, Ml, does not 
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show significant emission in the K band continuum, whereas the 
nucleus appears as the brightest peak in the K band continuum 
image but is not detected at MIR wavelengths. A very interest- 
ing peak is M7 to the Southeast, which is further discussed in 
section l4~4l 

3.2. SEDs of Clusters and Ring 

Figure [4] displays the various mid-IR spectra and photomet- 
ric measurements to construct the spectral energy distribution 
(SED) of the starburst ring of NGC7552. Here we compare 
our VISIR measurements with the literature data from Spitzer 
(IRAC, IRS) and ISO (SW S, PHOT). The ISOPHOT spectrum 
(ISiebe nmorge n et all 120041) was obtained through a 24" aper- 
ture, and the Spitzer-IRS spec trum from the SINGS Fifth Data 
Release dKennicutt et alJl2003l) integrated over a 27 "x43" wide 
spect ral map. In addition, ISO-SWS observations (IVerma et al.l 
[20031 not shown in Fig. |4]i revealed, within an aperture of 
14"x20", [Arm] and [Siv] line fluxes of 25 x 1(T 21 WcirT 2 
and 3 x 1(T 21 W cm , respectively. 

Within the given uncertainties, the agreement between the 
measured fluxes from different instruments and observing modes 
is very good. This may be surprising since the ISO and 
Spitzer measurements also include the nucleus of NGC7552 
whereas the ground-based measurements focus on the ring only. 
However, the contribution of the nucleus to the total emission 
from the central region becomes increasingly less toward longer 
wavelengths (cf. Figs. [3]b and [6]l and is negligible at mid-IR 
wavelengths. Furthermore, the fact that the Spitzer-IRAC pho- 
tometry agrees so well with the ISOPHOT measurements, and 
even the IRAS 12//m flux density (Fig.[U), suggests that essen- 
tially all of the mid-IR emission comes from the central 12" re- 
gion. For the total infrared luminosities of the individual clusters 
we refer to section POl 

3.3. Clusters versus Extended Emission 

Comparing the emission from the nine brightest mid-IR peaks 
to the total flux integrated over the ring, however, we find a 
significant mismatch. The large ratios between compact to total 
(compact plus diffuse) emission, quantified in Table [4] indicate 
that there must be substantial diffuse emission within the ring 
and outside the mid-IR peaks (XM,). The [Ne n] and 12/^m con- 
tinuum flux densities emitted by all identified clusters are only 
32% and 25%, respectively, of the emission integrated over the 
ring area. The contribution of the clusters to the total emission in 
the 8.6yum PAH band is even less, only 5%. Moreover, most of 
these 5% has its origin in the mid-IR source M7 (with less than 
1% contribution from the most luminous peak Ml). A substan- 
tial contribution of diffuse emission of [Ne n] has also been ob- 
served in the dwarf starburst galaxy NGC 5253, where ~ 80% of 
the [Ne n] line flux is diffuse and not directly associated with the 
central super star cluster (SSC) (iMartfn-Hernandezet al. 2005; 
Beir ao et al . 200&. 

We emphasise that ground-based observations in the mid- 
IR generally possess a much reduced sensitivity to low sur- 
face brightness. Here we use the measure d [Ne n] line fluxes 
for comparison: ISiebenmorgen et al.l(l2004l) derived (490 + 56) x 
10~ 21 WcirT 2 from TIMMI2 observations. IVerma et aT] (120031) 
measured 680 x 10~ 21 WcirT 2 from ISO-SWS data, and we de- 
rived (800±50)x 10~ 21 W cirT 2 for Spitzer-IRS from the 5" 1 data 
release of the SINGS Legacy Team. While the [Ne n] line fluxes 
from VISIR (Table [3]) and TIMMI2 agree reasonably well, the 



Table 4. MIR emission from clusters versus total emission in the 
ring. 





clusters / ring 


clusters / total 


ring / total 


[Nen] 


0.32° 


0.15 6 


0.48* 


12 fan cont. 


0.25° 


0.22* 


0.86 6 


8.6 PAH 


0.05" 


0.03 r 


0.65' 



"this work only; 6 i ncl. S INGS 5''' data release; r incl. data from 
ISiebenmorgen et al. (2004). 

space-based measurements are typically about a factor of two 
higher. In other words, the ground-based measurements under- 
estimate the emission from resolved low surface brightness re- 
gions. This technical shortcoming may be even more relevant for 
the derived PAH ratios. Hence, the above mentioned ratios from 
the first data column in Table |4] should only be considered an 
upper limit. 

In summary, a large fraction of the [Ne n] and PAH band 
emission is not directly associated with massive, young star clus- 
ters. In principle, the diffuse component may be due to an older 
population of stars (e.g., an aged generation of previous super 
star clusters), a more distributed mode of star formation (e.g., 
stars not localised in super star clusters), or a large fraction of 
radiation from the identified mid-IR peaks leaking out far into 
the surrounding interstellar medium. In any case, the localised 
emission diagnostics do not exclusively trace the most recent 
sites of star formation and need to be interpreted with care ( see 
sections I4~ni431 l. 

3.4. PAH Emission 

The spectra in Fig.|4]are dominated by strong PAH emission fea- 
tures at 6.2, 7.7, 8.6, 1 1.3 and 12.7 //m. PAHs are considered the 
most efficient specie s for stochastic, photo electric heating by UV 
photons in PDRs (Bakes & Tielenslll994l). They are commonly 
found in star form ing galaxies (e.g!7|Genzel & Cesarskv 2000; 
ISmifh et alJl2007h . However, intense UV fields ca n also lead to 
the g r adual destruction of PAH molecules (e.g ., iGeballe et al. 
1989; ICesarskvetaIlll99f3t iBeirao et al.l l2006h . PAHs may be 
considered as ove rall good tracers o f starburst activity in a sta- 
tistical sense (e.g. lBrandl et al.ll2006l and references therein). 

The low -resolution VISIR N band spectrum of peak Ml ob- 
tained is shown in Fig. [4] It is characterised by a rising contin- 
uum, typical for thermal emission by dust, a strong [Ne n] line, 
and a weak 11.3/zm PAH. We emphasise that the spectral flux 
density agrees very well with the photometric fluxes obtained 
through the VISIR filters, which agree well with the flux levels of 
the ISO and Spitzer spectra. Hence, we have good confidence in 
the absolute flux calibration of the VISIR spectra. Nevertheless, 
the weakness of the PAH emission feature in the spectrum is 
striking (see section l4~Tl i. 

Using PAHs as quantitative tracers of star formation is 
known to have some shortcomings. First, at high angular resolu- 
tion, the correspondence between observed PAH strength and the 
luminosity of the young clusters may break down, for reasons 
discussed in section |4~T1 Secondly, PAHs may also be excited 
in less UV-rich environments and can thus trace other sources 
besides massive you ng stars, such as planetary nebulae an d re- 
flection nebulae (e.g.ltlchida et alJl998tlLi & Dr aine 2002), and 
B stars dPeeters et al.l l2004b . Since these sources dominate the 
galactic stellar budget, i.e. are also much more numerous in the 
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Table 5. Continuum subtracted 8.6 and 11.3 jum PAH band 
fluxes in units of 1CT 21 WcirT 2 . 



Source 


8.6 PAH 


1 1.3 PAH 


Aperture 


Ml 


~0 


18" 


ff.' 7 


M2 


- 




a/7 


M7 


61 




a/ 7 


IMi 


~ 60 




a' 7 


Ring 


1130 




3"-8" 


TIMMI2 




330 


3" 


ISOPHOT 


1740 


1620 


24" 


IRS/Spitzer 




2570 


27"x43" 



Measurement obtained from the VISIR spectrum of Ml. 



volume that comprises the starburst ring in projection, we expect 
significant emission from these sources, in particular if the star- 
burst has already been going on for while and we do not observe 
the first generation of super star clusters any more. 

Table [5] lists the continuum subtracted PAH band fluxes 
for some mid-IR peaks and the integrated emission from the 
starburst ring. For comparison, measurements from TIMMI2, 
ISOPHOT and Spitzer-IRS are also included. 

3.5. Atomic Hydrogen Lines and Stellar Ages 

The information provided by the atomic hydrogen lines listed in 
Table [3] can be used in various ways. The intensity of the ex- 
tinction corrected Br-y line ( section \3.6\ can be used to calcu- 
late (see e.g., iHo et alj |l990) the number of hydrogen-ionising 
Lyman-continuum photons per second: 



A^Lyc = 8.7 x 10 4 



D 
kpc 



> Bry 



10~ 12 erg cm 2 s 



(1) 



From A^Lyc we can estimate an equivalent number of 07V 
stars (7Yq7 v). We use the cali bration of O star parameters pro- 
vided by Marti ns et al.l (|2005, Table 4) for the observational T e g 
scale and a number of log(7Yo7v) = 48.75 Lyman continuum 
photons per 07 V star. Both A^yc and Noiv are listed in Table [3] 

The equivalent widths of the hydrogen recombination lines 
Ha and Br-y are commonly used as age estimators of young stel- 
lar populations. These lines are predominantly produced by the 
most massive stars, which are short-lived. Hence, the strength of 
the line-to-continuum ratio is a strong function of cluster age as 
the most massive stars evolve. 

The evolutio n of the Br-y EW with time has been modelled 
in detail, e.g. bv iLeifherer et al. (1999, Starbust99). However, 
the relatively large distance to NGC 7552, the low EW of Br- 
y, and the dense ISM in the starburst ring suggest that modelling 
the MIR peaks as isolated young clusters may be insufficient. We 
have therefore included another step based on the combination 
of Starburst99 wit h the photo-ionisation code Mappings (e.g. 
Dopi taetal] 12000). which includes the radiation transfer in an 
evolving H n-region. iGroves et al.l (l2008al) have published a set 
of comprehensive, panchromatic starburst models, from where 
we took the equivalent width of Br-y as a function of the cluster 
age. We assume an instantaneous star formation, a thermal pres- 
sure log(P/k) - 5, and a so-called compactness logC = 5. This 
compactness parameter describes the time evolution of the H n- 
regi on, and is defined as logC = flog(^f t ) + |log(^y 
(see lGroves et all (l2008al) for details). 



To derive the cluster ages we compared the extinction cor- 
rected EWs of Br-y (TableO with the ages from the model. The 
resulting cluster ages are given in in the last column of Table [3] 

A couple of items are noteworthy considering our age es- 
timates. First, our approach is similar to the computed model 
tracks by Snii ders et al.l (120071 Fig. 14). Second, in comparison 
to "pure" Starburst99 estimates, we derive cluster ages that are 
systematically younger but only by 0.1 - 0.3 Myr. Third, the 
EW of Br-y is only accurate for young ages. Beyond 6.5 Myr, 
EW(Br-y) becomes very small, observational and model uncer- 
tainties begin to dominate, and for clusters older than 8 Myr Br- 
y is even observed in absorption and may reduce the emission 
features of nearby, younger clusters. Finally, the relative differ- 
ences among our derived cluster ages, however, can be consid- 
ered quite accurate: the difference between a cluster age of 5.5 
and 5.8 Myr corresponds to a difference in EWs of 18A and 8 A, 
respectively, which is much larger than the spectrophotometric 
errors. We shy away from quantifying formal errors on the in- 
dividual age estimates because we believe that the largest un- 
certainty may be given by the assumption that each photometric 
aperture includes only one, coeval stellar population. 

3.6. Extinction 

The N band includes the characteristic Si=0 stretching reso- 
nance of silicate-based dust, centred at 9.7 jum. However, since 
it is usually surrounded by strong PAH emission features to both 
sides (see Fig. |U left) its utilization as a quantitative measure 
of extinction in moderately dust enshrouded systems is very un- 
certain. On the other hand, the Br-y/Pa/S line ratio from the the 
SINFONI IFU data allows us to create an extinction map by 
comparison with the case B theoretic al ratio of 5.87 (assum - 
ing T e = 10 4 K and n e = 100 cm" 3 . iHummer & Storevll 19871) 
Assum ing a near-IR extinction law oc A~ l * dMartin & Whittetl 
1990) and R v = 3.1 dMatfaidl 19901) we calculated A v and list the 
values for the mid-IR peaks in Table [3] The average extinction 
of the nine peaks is Ay = 7.4 and reaches from 5.2 to values as 
high as 9.0 for source M3. Integrated over the entire ring area, 
the average extinction reaches "only" Ay = 5.3. 

With this information we constructed an extinction map, 
pixel-by-pixel for the central region of NGC 7552, which is 
shown in Fig. [8] These spatially resolved maps can be used to 
correct the observed Br-y line fluxes for extinction. The extinc- 
tion corrected Br-y line fluxes are listed in the third data column 
of Table [3] We note that there is, generally, a very good cor- 
respondence between the extinction peaks and our cluster posi- 
tions. 



4. Results and Discussion 

Generally speaking, there are four main reasons why MIR im- 
ages of starburst regions may reveal a different morphology than 
those taken at visible or near-IR wavelengths. All these effects 
should be kept in mind when discussing the properties of the 
"MIR clusters" to characterise the starburst ring. 

The first one is simply angular resolution which is propor- 
tional to oc A/D. Since the wavelength A is typically 5-10 times 
longer and the telescope diameter D is often smaller, at least for 
space telescopes, the resulting angular resolution of the MIR im- 
ages is much reduced. However, this effect does not play a role 
in this analysis - with the exception of the IRAC images pre- 
sented in Fig. [6]- based on MIR images of 073 - 0V4 resolution. 
The second one is simply the generally much lower sensitivity 
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Fig. 8. Extinction map computed from the Pa-/? and Br-y line 
emission maps. The mid-IR peaks from Fig.[3]a are indicated by 
red crosses. 



of ground-based MIR instruments. This effect will introduce a 
selection bias toward the most luminous clusters and overlook 
less luminous and older, more evolved clusters. The third one is 
dependency of the dust extinction on wavelength, which may re- 
veal heavily embedded clusters at MIR wavelengths while they 
are invisible at shorter wavelengths. The fourth effect is the fun- 
damental difference in the physical emission mechanism. While 
visible and NIR images mainly reveal the starlight of the cluster 
stars, the MIR emission is generally re-radiated light from the 
dust around the cluster. Fluctuations in the surrounding dust dis- 
tribution may thus affect the location of the MIR peak, and sub- 
structures of smaller clusters, surrounded by dust lanes may ap- 
pear as one large complex. In the following discussion we need 
to keep these issues in mind. 

4.1. Resolving the H n Region Structure 

In the most simplistic picture, the massive stars in a young clus- 
ter provide a sufficient number of photons with E y > 13.6 eV 
to ionise the surrounding hydrogen. The energetic radiation 
creates a so-called Hn region from where most of the emis- 
sion from the ionic states with higher excitation potential orig- 
inates: [Ne in] (40.96 eV), [Nen] (21.56eV), [S iv] (34.97 eV), 
[Sin] (23.34 eV), and [Arm] (27.63 eV). The Hn region is sur- 
rounded by molecular gas which is "only" exposed to far-UV 
radiation (6 - 13.6 eV ), which strongly influences its chemical 
and thermal structure (iTielens & Hoilenbachlll985l) . and which 
is responsible for most of the PAH emission. Beyond this photo- 
dissociation region (PDR) we expect the diffuse and partially 
neutral interstellar medium (ISM). In the most simplistic picture, 
the size of the H n region is given by the Stromgren radius 



Rs 



3 JVi 



Lyc 



[4nn 2 /3 B 



1/3 



(2) 



If we assume a density of 10 3 cm 3 , a recombination coef- 
ficient /3b = 2.6 x 10~ 13 cm 3 /s, and a Lyman continuum photon 



flux of N Lyc = 5 x 10 51 s" 1 (Table© we get R s = 5.4 pc or a 
region of about 10.8 pc in diameter in projection. 

In reality, however, the structure of giant PDR/H n regions is 
significantly more complex as ill ustrated by the closer ex amples 
of 30Doradus in the LMC (e.g..llndebetouw et al.ll2009l Fig. 1) 
and NGC 604 in M 33 (e.g. lHunter et al.lll996l Figs. 2 & 3). The 
main differences between the Stromgren picture and these gi- 
ant H ii regions are many-fold: they are not powered by a single, 
point-like cluster, but by several, distributed clusters; dust com- 
petes with the gas for ionising photons and will shrink the size of 
the H ii region; the radiation will not radiate perfectly isotropic 
and the interfaces between H n region, and PDR are also shaped 
by the stellar winds, outflows and supernovae. Hence, the above 
estimated Rs may only provide a rough estimate. In reality, the 
luminous PDRs are likely at larger distance from the centre of 
the H ii region. We would thus expect that spectra of the central 
cluster, taken at high angular resolution, will miss large contri- 
butions from the PDR. 

At the given distance to NGC 7552 of 19.5 Mpc and a slit 
width of 0775, the VISIR spectra cover a region of approxi- 
mately 70 pc in size. Indeed, the discrepancy in the strength 
of the 1 1.3 fim PAH feature between the high resolution VISIR 
spectrum (black line in Fig. HJ and the spatially integrated 
(and rescaled) Spitzer spectrum (red line in Fig. [4} is striking. 
Similarly, the ratios in Table|4]are based on photometric aperture 
sizes of 65 pc. While, within the same aperture size, the clusters 
contain 32% of the [Ne n] emission of the ring, they contain only 
5% of the PAH emission. In other words, zooming in on radial 
distances of 33 pc around the clusters reduces the observed PAH 
strength by a factor of six. This is a clear indication that we are 
resolving the H n/PDR complex. 

This may not come as a big surprise. The strong dependency 
of mid-IR spectral diagnostics on the aperture size has already 
been pointed out bv lMartfn-Hernandez et al.l (120061) . and a simi- 
lar effect was noted for the young super s tar clusters in t h e over - 
lap region of the A ntennae galaxies by ISniiders et al.l (120071) : 
Brandl et al . (2009). However, we want to emphasise two find- 
ings: First, a significant fraction of the [Ne n] emission escapes 
from the cluster regions. Second, at high angular resolution, the 
mid-IR spectra of the central clusters include only partial contri- 
butions from the PDR. 

Both findings have important general implications: On one 
hand, the physical properties of the young stellar populations 
(such as age or initial mass function, IMF) derived from lines 
fluxes integrated over the entire galaxy are very uncertain, as a 
large fraction of the line flux may not be associated with those 
populations. On the other hand, comparing spectral features that 
originate from physically and spatially different regions (such 
as ionic lines from H n regions and PAHs from PDRs) only for 
small apertures may lead to similar systematic errors. 

4.2. Cluster Luminosities and Star Formation Rates 

An important quantity to characterise starbursts is their bolomet- 
ric luminosity. In very dusty environments (Ay » 1) we assume 
that the bolometric luminosity Lboi is approximately given by the 
total (8 - 1000 /mi) i nfrared luminosity L ir. The latter is com- 
monly derived (e.g., ISanders et al.l 120031) from the four IRAS 
bands using the equation Lir = 562860 • D 2 ■ (13.485 n^m + 
5.\6S 25lim + 2.585 eotm + Sm^m), where S A is in Jy and D in 
Mpc. From the spatially integrated IRAS fluxes of NGC 7552 
in w e get L]r _ 8.5 x 10 10 Lq. Considering only the starburst 
ring, [Schinner er et al.l (119971) quote a bolometric luminosity of 
L bo i = 2.8 x 1O 1O L . 
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Of particular interest are the infrared luminosities of the indi- 
vidual clusters. However, observations that cover the FIR regime 
do not possess the required spatial resolution, and the VISIR ob- 
servations only cover 8-13 /urn. Hence, we bootstrap the VISIR 
fluxe densities measured through the 'NeII_2' filter (Table [TJ to 
Lir estimates of clusters from the literature. The 'NeII_2' fil- 
ter has been chosen as the longest wavelength filter of our data 
set, which d oes not not include significant emission or absorp- 
tion features. Brandl et al. 1 (120091 , Table 9) have provided SEDs 
and total infrared luminosities for six massive clusters in the 
Antennae galaxies. We have computed the ratios between the 
'NeIL2' filter band averaged and Lir for these six clus- 

ters. Adjusted for the distance of NGC 7552 we get the empirical 
relation 



Ljr = (30.1 ± 7.9) xF 12/Jffl , 



(3) 



where the 12 fim flux density is given in mJy and Lir in units of 
1O 6 L . The quoted uncertainty is the standard deviation of the 
six derived ratios. 

The estimated infrared luminosities are listed in Table [3] 
More than half of the total luminosity is provided by the clusters 
Ml, M2 and M7. Summing up the total luminosity of all MIR 
peak yields L^ = 2.1 x 1O 1O L . This is 75% of the bolomet- 
ric lu minosity of the starburst ring derived by ISchinnerer et alJ 
(fl997h . The good agreement of these two estimates, derived by 
completely different methods, supports our "bootstrapping" ap- 
pro ach. 

iKennicuttl fl998) has shown that the star formation rate 
(SFR) can be derived from Lir via: SFRtMoyr 1 ] = 4.5 x 
10~ 44 Lir [erg s ]. We note that this conversion strictly applies 
only to the continuous star formation approximation with ages of 
order 10100 Myr. Since our clusters are younger, the IR luminos- 
ity per uni t mas s of stars formed will be somewhat lower, and the 
iKennicutl d 1998h relation will overestimate the true SFRs. (This 
limitation would also apply to the common method of directly 
estimating the SFR from the IRAS fluxes). With the above con- 



\ ersion LjjJ s 2.1 x 1O 1O L corresponds to a SFR of 3.7M yr _1 



withi n only the main cluster s in the starburst ring. For compar- 
ison, ISchinnerer et al.l (1 19971) derived from the radio luminosity 
of the radio kno ts a star forma t ion rat e of 1M yr~' in the ring of 
NGC 7552, and lCalzetti et al.l d2010l) derived from Spitzer pho- 
tometry taken by the SINGS team lOJMpyr' 1 for the entire 
galaxy . The sample of 22 nuclear ring galaxies of Mazzuca et al. 
(2008) covered SFRs from O.1M yr _1 to 9.9M yr -1 , with a me- 
dian value of 2.2M yr~' . Our value of 3.7M yr~' fits very well 
within these numbers. 

4.3. Cluster Ages 

The ages of clusters Ml - M9 are listed in the last column of 
Tableland lie between 5.5 and 6.3 Myr with a mean age of 
5.9 + 0.3 Myr. Although they have been derived under the sim- 
plifying assumption of instantaneous star formation we can draw 
two important conclusions: First, there are no luminous young 
clusters associated with the mid-IR peaks that indicate very re- 
cent (<k 5 Myr) star formation. Second, the age spread between 
the clusters associated with the mid-IR peaks is relatively small. 

To investigate whether high extinction may hide the 
youngest clusters even at K band, we also use independent di- 
agnostics at MIR wavelengths. A commonly used age diagnos- 
tic is the ratio of the two neon fine structure lines, [Ne m] and 
[Nen], because it is largely independent of density and extinc- 
tion effects. With excitation potentials of 40.96 eV and 21.56 eV, 



respectively, the ratio of [Ne m] / [Ne n] measures, to first degree, 
the hardness of the radiation field, which is mainly dominated by 
the massive O stars and thus a strong function of age and initial 
stellar mass function (IMF) of an instantaneous starburst. (We 
note that factors other than age may also play a role in affecting 
the observed line ratios, such as high metallicity ( via stellar evo- 
lution and line blanketing), high electron density (Sniider s et al.l 
2007), and variations of the upper IMF). Unfortunately, we have 
no measure of the [Nem]15.56//m line, which lies outside the 
atmospheric transmission window, but we can still estimate the 
MIR fine structure line ratios in two different ways: 

First, [Ne in] has been measured from space. Both ISO-S WS 
and Spitzer-IRS found rather low [Ne in] / [Ne n] « 0.08 through 
large apertures that covered the central region. The large dis- 
crepancy between compact and extended emission, as discussed 
in section 1331 complicates the interpretation of the line ratios. 
However, if we assume that most of the [Ne in] is only produced 
by the massive O stars in the clusters - which is not an unre- 
alistic assumption given the high excitation potential of [Ne in] 
— we can compute th e ratio between the [Ne in] flux given by 
iThornlev et al.l d2000l) and the sum of [Ne n] from all clusters 
(Table 0. Under this assu mption, we derive [N e ml / [Ne nj « 
0.41. Starburst models (e.g. Sniid ers et al.l l2007) indicate that a 
ratio of 0.4 is in good agreement with an age of approximately 
5 Myr. 

Second, for the brightest MIR peak Ml we obtained a VISIR 
spectrum that covers also the [S iv] 10.5 1 fim line. Unfortunately, 
the [Siv] line was not detected in the VISIR spectrum of Ml, 
so we need to follow a different route. [S iv] has an excitation 
potential of 34.97 eV, which is cl ose to the 40.96 eV of [Nem], 
and may thus serve as a substitute. iGroves et al.l(l2008bl) provide 
an empirical calibration of [S iv] / [Ne n] against [Ne in] / [Ne n], 
parametrised as: 



log 



[Ne in] 
[Neil] 



a log 



[Siv] 
[Ne II] 



(4) 



From fits to a large s ample of 97 extragala ctic H n regions and 
56 starburst galaxies IGroves et al.l (1200 8b) derived the param- 
eters a = 0.82 (0.65) and (3 = 0.36 (0.32) for extragalactic 
Hn regions (starburst galaxies). From the [Nen] flux listed in 
Table [3] the above correlation(s), and [Ne m] / [Ne n] a 0.41, 
we would expect a [Siv] line flux from peak Ml of L^ 1 ^ = 

3.1 x 10~ 21 Wcnr 2 (Ff** = 2.1 x 10~ 21 Wcnr 2 ). In both cases, 

lo 1VI 

the line is too weak to be detected in our VISIR observations. Of 
course, this estimate assumes that the properties of Ml are sim- 
ilar to those of the other clusters in the starburst ring. 

Independently, we can derive an upper limit on the [S iv] line 
flux from the noise in the 10.4 - 10.6 //m range. This yields a 3cr 
upper limit of L^ v] = 2.3 x 10~ 20 WcirT 2 . Hence, we derive 

an upper limit on the ratio of f ffi V ]/f[Neii ] ^ 0.12. Comparing 
this ratio with the models of Snii ders et al.l (120071 Fig. 9) yields 
a lower limit of 4 Myr on the age of peak Ml, consistent with 
the Br-y age estimate. 



4.4. The Intriguing Peak M7 

A rather unusual and thus very interesting source is peak M7, 
located to the Southeast of the center. A qualitative comparison 
between the three prominent peaks, Ml, M7, and the nucleus, 
is given below, where '+' means bright and '-' means faint or 
undetected: 
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Kcont 8 .6 /mi PAH [Nen] 

Nucleus + 

Ml - - + 

M7 + + + 

M7 is equally prominent in Br-y, [Ne n], and the K band con- 
tinuum. What is most remarkable, however, is its strong 8.6 /mi 
PAH emission. Fig [2] (center) shows that M7 is the brightest 
8.6/im PAH emitter in the center of NGC7552. Furthermore, 
M7 is a peculiar source in two more regards. 

First, the velocity field of the nuclear region (see section l4~8l 
and Fig. O, while overall relatively symmetric, reveals a distinct 
"nose" in the contour line of +50kms~l, approximately 3" to 
the Southwest, where M7 is located. Apparently, the velocity 
field in the r egion of M7 is sl i ghtly d isturbed. 

Second, Rosenberg et al.l (1201 2l) have recently investigated 
the relationship between the [Fe n] 1 .26/mi luminosity and the 
supernova rate in a sample of 11 nearby starburst galaxies, 
including NGC7552. They derived the supernova rate from 
Starburst99 modelling of the Br-y emission on a pixel-by-pixel 
basis and found a very good correlation across all 1 1 galaxies be- 
tween the derived supernova rates and the measured [Fe n] line 
luminosities - with one exception, namely M7 in NGC7552. 
According to lRosenberg et al.l(l2012l their Fig. 7), the supernova 
rate derived from the Br-y emission exceeds by far the corre- 
sponding [Fe n] intensity. 

We note that the derived age of M7, its dust extinction, Br-y 
equivalent width, and total luminosity are well within the range 
of parameters of the other peaks (Table |3}- Nevertheless, peak 
M7 is not likely a simple, single massive cluster (see section l4~6l 
for a discussion of possible causes.) 

4.5. Comparison with other Massive Clusters 

The main question here is: Are the MIR selected clusters in 
NGC 7552 in any way special, or similar to massive young clus- 
ters in other starburst galaxies? This can be best addressed by 
comparing the cluster properties in Table [3]with several "refer- 
ence clusters" from the literature (Table |6). 

While the globally integrated ratio of [Ne m] / [Ne n] « 0.08 
is significantly below the mean value of 20 starburst galaxies, 
the "corrected" value of 0.41 (section 03) places NGC 7552 
well within the range of starbursts and only a factor of two be- 
low the mean radiation hardness of Galactic H n regions. The 
real outlier in this regard in Table |6]is NGC 3603, a very young 
Galactic Hn region, dominated by a single compact cluster, 
and NGC 5253. For the latter it was shown that the bulk of the 
high-excitation [S iv] and [Ne m] fine structure line emission is 
also associated with the single compact cluster C2, whereas the 
diffuse component - showing [ Arn~|, fSml, and [Nenl emis- 
sion - is much m ore extended (Mart in-Hernandez et al.l 12005; 
Beirao et al. 2006). In terms of Lyman continuum photons A^yc, 
our peak Ml is only a factor th ree below the radio super - nebul a 
NGC 5223 C2 (2 x 10 52 s _1 ; iMartm-Hernandez et all l200l . 
However, in terms of the ratio between ionising photons and 
the surrounding PDR mass, the peaks in NGC 7552 resemble 
the massive clusters in more luminous starburst galaxies, while 
NGC 3603 and NGC 5253 may be better compared to blue com- 
pact dwarf galaxies. 

The total infrared luminosities of the peaks Ml through 
M9 are, on average, somewhat fainter than the most lumi- 
nous Antennae cluster, but agree well with e.g., NGC 5223 C2, 
NGC 1365 M6 and NGC 4038/39 peak 5, which display simi- 
larly high extinction. 



Since we have no good direct tracer of stellar mass (e.g. 
dynamical masses) we refrain from providing cluster mass es- 
timates here. However, a comparison with the reference clusters 
of similar ages, Lir and A^yc in Table [6] suggests that our clus- 
ters have typical masses of a few million solar masses, which 
is not extraordinarily large. Significantly more massive young 
clusters, up to 3 _xl0 1 M G , have been found e.g. in NGC 1365 by 
iGalliano et all d2008l) . 

The ages of the clusters in NGC 7552 are also well within 
the range of the MIR luminous clusters in Table [6] It is evident 
that most of the massive clusters listed here managed to retain 
a substantial amount of extinction after 5-7 Myr. An excep- 
tion here is Arp 143. We note, however, that the starburst ring in 
Arp 143 results from a shock wave triggered by the recent head- 
on co llision between the two galaxies, NGC 2444 and NGC 2445 
(e.g jBeirao et al.l2009l and references therein). This scenario is 
physically different from dynamical instabilities in a quasi-stable 
nuclear region, and this difference will be reflected in the den- 
sity of the ISM and its gas supplies (see also section 14.71) . In 
summary, we conclude that the clusters in the ring of NGC 7552 
appear like typical unit cells of starbursts. 

4.6. "Clusters" and "Peaks" may be Misleading Terms 

So far we have, in a very simplistic way, equated a MIR "peak" 
with a super star cluster. Both species are closely related in star- 
bursts, but there may not necessarily be a direct correspondence. 
In particular, one MIR peak may consist of several smaller clus- 
ters. We remind the reader that our angular resolution at near- 
and mid-IR wavelengths corresponds to about 30 pc, which is 
about one to two orders of magnitude lar ger than the typica l 
core radius of a massive cluster (e.g., R136. lBrandl et al.l fl996). 
On one hand, it is likely that the hierarchical cloud fragmenta- 
tion in a dense and turbu lent ISM will lead to a conglomerate of 
smaller clusters (e.g.. iBonnell et aLll2003l) . On the other hand, 
once a massive cluster has formed in a dense environment it 
may trigger the f ormation of other clusters in its vicinity (e.g., 
Deharveng et al. 2005). While both scenarios impact dynamical 
mass estimates and cluster evolution, the latter scenario has the 
most severe effect as the peak cannot be any longer considered a 
coeval population, formed in an instantaneous event. 

Observationally, there i s strong evidence fo r a more com- 
plex scenario. For instance, iBastian et al 1 (120091) found that the 
youngest clusters in the Antennae are not isolated but part of an 
extended hierarchy of star-forming regions. In those "nuclei" of 
star formation, one cluster is often surrounded by other clusters, 
which cannot be resolved with the current combination of tele- 
scopes an d wavelengths. In th e overlap region of the Antennae 
galaxies, Sniiders et al. (2007) found on scales below the reso- 
lution limit a complex structure, which contains several young 
stellar clusters embedded in clum py gas. Simi l arly, i n a statisti- 
cal analysis of Antennae clusters Meng el et al.l(l2005l) found that 
the 10 Myr old clusters with high extinction are preferentially lo- 
cated near younger clusters. This suggests that in environments 
where the gas is not rapidly removed, ongoing - maybe even 
triggered - cluster formation is likely to happen. 

These aspects are very relevant to the derived cluster proper- 
ties as an unresolved complex of (sub-)clusters with significant 
age gradient will appear differently. While the high excitation 
lines (e.g. [Nem]) will be dominated by the young cluster, the 
lower excitation lines (e.g. [Ne n]) will primarily originate in the 
older components, lowering the combined neon line ratio. In ad- 
dition, the IMF - when normalised to the most luminous stars 
- will appear steeper (bottom heavy). Furthermore, dynamical 
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Table 6. Data of selected luminous clusters from the literature for comparison 





Galactic 


Starburst 


NGC 3603 c 


NGC5253 rf 


NGC 1365 e 


NGC4038/39 f 


Arp 143 s 


M82 A 




H ii regions" 


galaxies* 




C2 


M5 


M6 


peak 1 


peak 5 


knot E 


knot F 


region 2 region 3 


L m [10"L o ] 




50 


0.01 


1.0 




0.8 


11.4 


3.9 








M [10 6 M o ] 








» 1 


38 


25 


3 




25 


6 




A v [mag] 




4.3 


< 1 


17 


3.2 


8.5 


4.5 


10.5 


1.2 


1.3 


^15 * 25 


Age [Myr] 






~ 1 


< 6 


7 


7 


2.5 


4.5 


3.5 


2.5 


5.5 5.5 


[Nem] / [Nen] 


0.8 


0.29' 


3.13 


8.4 






0.63 


0.69 


0.52 


0.67 


0.13 0.19 



' Giveon et af] d2002h iThornlev et al.l (120001). 
( 2008) ^Brandl et alj f2009) * Beir ao et al.U2009nlBeirao et al 



Brandl et al.l d2006l) Cebouteiller etail d2008l) iMartfn-He mandez et al. (2005) 'Gallian oetail 
<2008h ' mean of 20 starburst galaxies, excluding the dwarfs II Zw 40 and NGC 5253. 



mass estimates yield incorrect results since all components of 
the cluster complex will contribute to the diagnostic line but are 
not in virial equilibrium. And finally, significant patches of dust 
may be trapped in lanes between sub-clusters and remain asso- 
ciated with the observed peak for a longer time, as discussed in 
the next section. 



4. 7. Cluster Evolution and Extinction 

Table [3] shows that the clusters of NGC 7552 have an average 
extinction of 7.4 mag (at an average age of 5.9 Myr). The large 
amount of gas and dust, still associated with our MIR peaks af- 
ter relatively long time, indicates that the gas and dust disper- 
sion has been rather inefficient during the first six million years. 
However, selecting clusters from MIR data introduces a selec- 
tion bias toward the most heavily extincted clusters. This can be 
best illustrated for the Antennae galaxies: 

The optically selected clusters within the Antennae have 
extinctions of < 4 mag (with a mean o f 2.6 mag in the age 
range 0-4 Myr (Whitmo re & Zhang|2002l). In a study of K band 
selected clusters in the Antennae (Menge let alj l200"5h found 
strongly variable extinction with an average value of Ay = 
1.3 mag. However, the MIR-selected clusters in the dusty over- 
lap region of the Antennae have ex tinctions of Ay > 4 up to 
Ay ~ 10 mag dBrandl et al.l 120091) . Although the latter may 
not dominanate by num bers, they are no exotic species either: 
IWhitmore etail d2010l) found that 16% + 6% of the IR-bright 
clusters in the Ant ennae are stil l heavi ly obscured with values 
of Ay > 3 mag. iMengel et al.l ((2005) found a clear trend of 
lower Ay with increasing age . Several studies (e.g. Mengel et al. 
120051: IWhitmore et al.ll20 10) suggest that the typical time scale 
for massive clusters to emerge from their natal dust cocoons is 
less than 10 Myr. 

However, the local environment and the location of the clus- 
ters within the galaxy is likely to play an important role in the ef- 
ficiency of the dust removal. It is conceivable that the expansion 
of the H n region goes faster in lower density environments (e.g., 
in dwarf galaxies or the outer spiral arms) than in higher density 
regions of galactic nuclei where the gas and dust are trap ped in 
a pote ntial and even get continuously replenished. Mengel et al. 
(2005) argue that an older cluster which suffers high extinction 
may no longer be directly surrounded and obscured by its own 
dust cocoon, but just happens to be located in a denser region of 
more recent star formation with higher dust content. 

While this picture is likely to apply to clusters in star- 
burst regions in general, and to clusters in the starburst ring of 
NGC 7552 in particular, we would like to reemphasize the "spe- 
cial" and distinct nature of MIR-selected "peaks". Several stud- 



ies (e.g., IWhitmore & Zhang||2002t IWhitmore et al.ll2010l) have 
found that, by number, the vast majority of hundreds of clus- 
ters in the Antennae is optically visible and can be well stud- 
ied with the HST On the other hand, the few highly reddened 
MIR peaks in the overlap region, although almost insignificant 
by number counts in comparison, account for approximately half 
of the total lumino sity - and thus o f the total star formation - 
of the Antennae (Bran dl et aDl2009l) . It is likely that the latter 
species does not exist in galactic disks and can only be found in 
the densest environments. However, the study of clusters in star- 
burst nuclei, ulta-luminous and sub-millimeter galaxies requires 
infrared observations at high angular resolution (cf. Fig. |7J. 

4.8. Star Formation History in the Ring 

In the Introduction (section [TJ we have summarised the contro- 
versy around the physical origin and location of starburst rings. 
To complicate matters, the starburst ring does not appear to be an 
narrow ring with a series of young clusters lined up on a circle. 
Instead, we observe (e.g., Fig. [3} an annular region of ~ 200pc 
width with increased and non-uniform star formation activity. 

Here we take the existence of a ring for granted and want to 
investigate how and where the clusters form within that ring. 
We assume that the gas gets efficiently trans ported from the 
outer regions to th e centre via the nuclear bar (ISchwarzlll98"Tl 
[Ellison et al. I l201ll) . The gas moves radially inward along the 
bar dust lanes and accumulates at the so-called contact points, 
the intersections between the dust lanes and the ring. These 
contact points are expected to be approximatel y perpendicular 
to the position angle of the bar major axis dRegan & Teubenl 
120031) . Since the orientation of the major axis of the bar is 
roughl y East- West while the c ontact points are to the North and 
South dSchinnerer et al.l ll997). the observations support the the- 
ory. However, the details of where, how, and when the gas gets 
turned into stellar clusters are still uncertain and could possibly 
be explained with one of the following two scenarios. 

In the first scenario, gas accumulates gradually in the nuclear 
region along the ring. Once a critical density is reached locally, 
the gas becomes unstable to gravitational collapse, and shocks 
may t rigger the formation of a massive star cluster (Elmegreen 
1994) or lead to further cloud fragmentation due t o supersonic 
turbulence (e.g jKlessenl l2001l;lBonnell et al.| [2003), resulting in 
an ensemble of clusters at that location. Since these instabilities 
occur randomly with respect to the galactic structure, no system- 
atic age gradients or patterns within the ring would exist. 

In the second scenario, most of the gas enters the ring at the 
intersection between the bar and the ring. Shocks play an im- 
portant role here to get rid of the angular momentum, which 
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is essential to mov e the gas into the central ring (see e.g., 
Mazz uca et ai] l2008l and references therein). The pile-up of gas 
at these two contact points can be expected to ignite star for- 
mation in short-lived bursts. The clusters formed here will orbit 
along the ring, and age as they travel along the ring. In the mean- 
time, new starburst events near the contact points may produce 
new generations of young star clusters. Hence, observationally, 
one would expect to find the youngest clusters of the ring near 
the contact points and increasingly older clusters in the direction 
of rotation. 

Some e v idence for the latter scenario was provided by 
iBoker et aD (|2008) for several galaxies, and in particular for 
NGC 613. In a more comp rehensive study of 22 barred galax- 
ies , Ma zzuca et al. (2008) found that 10 out of 22 ring galax- 
ies showed random variations in their H n region properties with 
no apparent age gradient. However, 9 (out of 22) of their rings 
showed partial gradients along 25% - 30% of their rings. Only 
three galaxies (NGC 1343, NGC 1530 and NGC 4321) showed a 
clear bipolar age gradient along the ring. iMazzuca et al .1 (120081) 
could not relate the presence or absence of age gradients with 
the morphology of the rings or their host galaxies. 
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Fig. 9. Velocity field in [km/s] of the nuclear region of NGC 
7552 derived from the 2. 166 fim Br-y line shifts. 



In order to investigate a possible age gradient across the 
clusters in NGC 7552 we constructed a velocity map from 
the SINFONI Br-y data (Fig. [9]). The fine offsets were con- 
verted into rotational velocities f or a recessional v elocity of 
NGC7552 of (1640 + 12)kms- 1 dMould et al.ll2000l including 
"Virgo+GA+Shaply", ). With this velocity offset the kinematic 
centre position appears to lie close to the nucleus of NGC 7552 
as traced by the K band continuum peak, and the velocity field 
appears rather symmetric to both sides of the centre. The zero 
velocity line of the Br-y image is at the position angle of ~ 9° 
with the kinematic major axis at ~ 99°. The physical velocity of 
a ring element is then related to the observed, Doppler-shifted 
velocity via: 

_ Vobs 

sin(i) cos(#) 



where i is the inclination of 28° dFeinstein et al.l fl990l) and 6 
is the angle in the plane of the ring measured from the posi- 
tion angle. From Fig. [9] we derive a mean v r j ng « 150km s _1 , 
which is in perfect agreement with t he typ ical rotation velocity 
of v ro t = 150 km s _1 IMazzuca et alj d2008) found for their sam- 
ple of nuclear rings. 

We can now estimate the travel time of a cluster at r — 275 pc 
along a quarter (90°) segment of the ring via: 

Wi = ^ = !^^*2.8Myr (6) 
Vring 150km /s 



| Contact / 
; point J 




| point « 



Fig. 10. Schematic outline of the cluster age distribution with 
respect to the location of the contact points in NGC 7552. The 
green circles indicate the positions of the MIR peaks and the 
inscribed numbers their ages in million years. The direction and 
time scale of rotation are also indicated. 

Fig.[10]shows a simplified representation of the scenario for 
NGC 7552 with both cluster locations and ages, and the travel 
timescales indicated. The youngest knots Ml and M4 appear in- 
deed to be located at or near the Northern contact point, and the 
oldest knots M6, M8 and M9 are located more than 90° away 
from the contact points. This signature provides evidence that 
the cluster age distribution is not completely random but some- 
how connected to the co ntact points of the ring . Our finding also 
agrees with the result of Ma zzuca et al. (2008) that in two-thirds 
of their 22 barred galaxies the location of the youngest H n re- 
gion is within 20° of the contact points. 

Despite the limited number of MIR peaks, we find some 
weak evidence for an age gradient in the direction of the ring 
rotation. However, the gradual age difference of the clusters is 
much smaller than the estimated cluster travel time. The pres- 
ence of older clusters and the lack of very young clusters close 
to the contact points raise some doubts in this simple (second) 
scenario. Our observations clearly do not support the hypothe- 
sis that the high gas mass inflow rate will immediately trigger 
gravitational collapse and ignite massive cluster formation. 
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However, since the cluster ages are relatively large, they are 
likely to have travelled a significant distance from their birth- 
place to their observed location. It is intriguing that the age of 
peak Ml, now located near the northern contact point, is exactly 
the travel time from the southern contact point to its current lo- 
cation (half a ring orbit). It appears therefore possible that a dra- 
matic episode of gas inflow at the contact points about 5.6 Myr 
ago has triggered massive cluster formation at those locations. 
However, while this scenario would explain the location of at 
least some of the MIR peaks, it would also create a new prob- 
lem: why should the gas inflow have stopped 5.6 Myr ago, at 
least to a level where no more massive young clusters are being 
formed? Since we have no additional evidence for such a sce- 
nario we don't consider it very likely. 

In summary, the central cluster formation in NGC7552 is 
not dominated by a simple bar/ring geometry, but it is not com- 
pletely decoupled from the ring mo tion either. Eithe r some neg- 
ative feedback effects are at play dDale et al .1120051) or the gas 
density will be built up more gradually over a larger ring seg- 
ment, and the local physical conditions become dominant. 

5. Summary 

We have observed the starburst ring galaxy NGC 7552 with two 
instruments on ESO's VLT, namely the Imager and Spectrograph 
for the mid-IR (VISIR) in both imaging and spectroscopy mode, 
and the Spectrograph for INtegral Field Observations in the Near 
Infrared (SINFONI) at the VLT. The angular resolution of the 
VISIR observations is 073 - 074, very close to the diffraction 
limit of VLT at 10/mi. We compared these data to ISO and 
Spitzer observations. Within the given uncertainties, the agree- 
ment between the measured fluxes from the different instruments 
and observing modes is very good. Our main findings can be 
summarized as follows: 

- The starburst ring is clearly identifiable at MIR wavelengths. 
The [Nen]12.8jum line emission correlates well with the 
NIR Br-y emssion. 

- We find an excellent correspondence between the [Nen] 
emission regions and the highly dust extincted dark ring in 
optical HST observations. The latter look very similar to the 
K band continuum. 

- We have identified nine peaks of unresolved MIR emission. 
Our strongest 12 /mi source does not show significant emis- 
sion in the K band continuum. 

- We do not detect MIR emission from the nucleus of 
NGC 7552, which is very prominent at optical and NIR con- 
tinuum wavelengths. 

- The multi-wavelength data indicate that the commonly used 
term 'starburst ring' is insufficient to properly characterize a 
rather complex picture of cluster formation within an annular 
region of more than 100 parsec width. 

- A large fraction of the [Ne n] and PAH band emission is not 
directly associated with the location of the massive, young 
star clusters. At the cluster positions we measured about 32% 
of the total [Nen] emission of the ring but only 5% of the 
PAH emission. 

- A comparison of the strength of the 11.3/im PAH feature 
between VISIR and Spitzer spectra indicates that zooming 
in on the clusters on scales of 33 pc reduces the observed 
PAH strength by a factor of six. This can be explained by 
spatially resolving the H n region from the PDR complex. 

- The individual cluster properties, in terms of infrared lumi- 
nosities, ages, Lyman continuum photons, and extinction, are 



in good agreement with the properties of other MIR-selected 
massive clusters in other galaxies. 

- From the ratio of Br-y/PayS we computed an extinction map. 
The average extinction of the nine peaks is Ay =7.4 and 
reaches from 5.2 to 9.0. While the average extinction of our 
MIR-selected clusters is more than five times higher than that 
of K band selected clusters in the Antennae, it is similar to 
MIR-selected clusters in other dense starburst environments. 
Apparently, the mechanism of gas and dust removal is less 
efficent in these extreme environments. 

- The total luminosity of the nine MIR peak is L/r = 2.1 x 
10 10 L o , which is about 75% of the bolometric luminosity of 
the starburst ring. In the continuous star formation approxi- 
mation this would correspond to a SFR of 3.7M yr _1 . 

- We determined the cluster ages from the equivalent width 
of the Br-y line. They lie within the range of 5.9 + 0.3 Myr. 
Independently, we estimated a ratio of [Ne m] / [Ne n] « 0.4 1 , 
which is in good agreement with an age of approximately 
5 Myr. 

- The youngest massive clusters are located near the Northern 
contact point of the ring while the oldeer peaks are ob- 
served further away from the contact points. However, the 
age spread among the clusters of 0.8 Myr is small compared 
to the travel time of ~ 5.6 Myr for half an orbit within 
the starburst ring. We cannot rule out that the clusters were 
formed near the contact points more than 5 Myr ago, but we 
find no strong evidence for the scenario where the continu- 
ous inflow of gas leads to the ongoing formation of massive 
clusters near the contact points. It appears more likely that 
the gas density build up more gradually over larger ring seg- 
ments, and that the local physical conditions then determine 
the cluster formation. 

- Even at the resolution of VISIR on an eight-meter VLT, we 
cannot rule out that an unresolved MIR peak is not a sin- 
gle super star cluster but may consist of several sub-clusters, 
possibly with some age spread. If so, this would affect the de- 
rived ages, the slope of the IMF, and dynamical cluster mass 
estimates. This issue can only be solved with milli-arcsecond 
resolution at MIR wavelengths as will be provided by the 
next generation of extremely large telescopes (ELTs). 

Acknowledgements. We would like to thank Brent Groves for helpful discus- 
sions concerning the cluster age estimates, and Brad Whitmore for a discussion 
on the observational aspects of dust extinguished clusters. We would also like to 
thank the anonymous referee for spotting an error in the NIR data reduction in 
the draft version. 



References 

Bakes, E. L. O. & Tielens, A. G. G. M. 1994, ApJ, 427, 822 
Bastian, N., Trancho, G., Konstantopoulos, I. S., & Miller, B. W. 2009, ApJ, 701, 
607 

Beck, R. 2011, ArXiv e-prints 

Beck, S. C., Turner, J. L., Langland-Shula, L. E., et al. 2002, AJ, 124, 2516 
Beirao, P., Appleton, P. N., Brandl, B. R., et al. 2009, ApJ, 693, 1650 
Beirao, P., Brandl, B. R., Appleton, P. N. andGroves, B., et al. 2008, ApJ, 676, 
304 

Beirao, P., Brandl, B. R., Devost, D., et al. 2006, ApJ Lett., 643, LI 

Binney, J. & Tremaine, S. 2008, Galactic Dynamics: Second Edition, ed. Binney, 

J. & Tremaine, S. (Princeton University Press) 
Boker, T., Falcon-Barroso, J., Schinnerer, E., Knapen, J. H., & Ryder, S. 2008, 

AJ, 135, 479 

Bonnell, I. A., Bate, M. R., & Vine, S. G. 2003, MNRAS, 343, 413 
Bonnet, H., Abuter, R., Baker, A., et al. 2004, The Messenger, 1 17, 17 
Brandl, B., Sams, B. J., Bertoldi, R, et al. 1996, ApJ, 466, 254 
Brandl, B. R., Bernard-Salas, J., Spoon, H. W. W., et al. 2006, ApJ, 653, 1129 
Brandl, B. R., Snijders, L., den Brok, M., et al. 2009, ApJ, 699, 1982 
Buta, R. & Combes, R 1996, Fund. Cosmic Phys., 17, 95 



15 



B. R. Brandl et al.: The Starburst Ring in NGC7552 



Calzetti, D., Wu, S.-Y., Hong, S., et al. 2010, ApJ, 714, 1256 

Cesarsky, D., Lequeux, J., Abergel, A., et al. 1996, A&A, 315, 309 

Dale, D. A., Gil de Paz, A., Gordon, K. D., et al. 2007, ApJ, 655, 863 

Dale, J. E., Bonnell, I. A., Clarke, C. J., & Bate, M. R. 2005, MNRAS, 358, 291 

de Vaucouleurs, G, de Vaucouleurs, A., Corwin, Jr., H. G, et al. 1991, 

Third Reference Catalogue of Bright Galaxies (Volume 1-3, XII, 2069 

pp. Springer- Verlag Berlin Heidelberg New York) 
Deharveng, L., Zavagno, A., & Caplan, J. 2005, A&A, 433, 565 
Dopita, M. A., Kewley, L. J., Heisler, C. A., & Sutherland, R. S. 2000, ApJ, 542, 

224 

Durret, F. & Bergeron, J. 1988, A&AS, 75, 273 

Ellison, S. L., Nair, R, Patton, D. R., et al. 201 1, ArXiv e-prints 

Elmegreen, B. G. 1994, ApJ, 425, 73 

Feinstein, C, Mendez, M., Vega, I., & Forte, J. C. 1990, A&A, 239, 90 
Forbes, D. A., Kotilainen, J. K., & Moorwood, A. F. M. 1994a, ApJ Lett., 433, 
L13 

Forbes, D. A., Norris, R. P., Williger, G. M., & Smith, R. C. 1994b, AJ, 107, 984 

Galliano, E., Alloin, D., Pantin, E., et al. 2008, A&A, 492, 3 

Geballe, T. R., Tielens, A. G. G. M., Allamandola, L. J., Moorhouse, A., & 

Brand, P. W. J. L. 1989, ApJ, 341, 278 
Genzel, R. & Cesarsky, C. J. 2000, ARA&A, 38, 761 

Gil de Paz, A., Boissier, S., Madore, B. F., et al. 2006, ArXiv Astrophysics e- 
prints 

Gillessen, S., Davies, R., Kissler-Patig, M., et al. 2005, The Messenger, 120, 26 
Giveon, U., Sternberg, A., Lutz, D., Feuchtgruber, H, & Pauldrach, A. W. A. 

2002, ApJ, 566, 880 

Gorjian, V., Turner, J. L., & Beck, S. C. 2001, ApJ Lett., 554, L29 

Groves, B., Dopita, M. A., Sutherland, R. S., et al. 2008a, 176, 438 

Groves, B., Nefs, B., & Brandl, B. 2008b, 391, LI 13 

Ho, P. T. P., Beck, S. C, & Turner, J. L. 1990, ApJ, 349, 57 

Hummer, D. G. & Storey, P. J. 1987, MNRAS, 224, 801 

Hunter, D. A., Baum, W. A., O'Neil, E. J., & Lynds, R. 1996, ApJ, 456, 174 

Indebetouw, R., de Messieres, G. E., Madden, S., et al. 2009, ApJ, 694, 84 

Kenney, J. D. P., Carlstrom, J. E., & Young, J. S. 1993, ApJ, 418, 687 

Kennicutt, Jr., R. C. 1998, in Astronomical Society of the Pacific Conference 

Series, Vol. 142, The Stellar Initial Mass Function (38th Herstmonceux 

Conference), ed. G. Gilmore & D. Howell, 1 — h- 
Kennicutt, Jr., R. C, Armus, L., Bendo, G, et al. 2003, PASP, 1 15, 928 
Klessen, R. S. 2001, ApJ, 556, 837 
Knapen, J. H. 2005, A&A, 429, 141 

Lagage, P. O., Pel, J. W., Authier, M., et al. 2004, The Messenger, 117, 12 
Lebouteiller, V, Bernard-Salas, J., Brandl, B., et al. 2008, ApJ, 680, 398 
Leitherer, C, Schaerer, D., Goldader, J. D., et al. 1999, ApJS, 123, 3 
Li, A. & Draine, B. T. 2002, ApJ, 572, 232 
Liu, J.-F. & Bregman, J. N. 2005, ApJS, 157, 59 
Martin, P. G. & Whittet, D. C. B. 1990, ApJ, 357, 113 

Martin-Hernandez, N. L., Schaerer, D., Peeters, E., Tielens, A. G. G. M., & 

Sauvage, M. 2006, A&A, 455, 853 
Martin-Hernandez, N. L., Schaerer, D., & Sauvage, M. 2005, A&A, 429, 449 
Martins, F, Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 
Mathis, J. S. 1990, ARA&A, 28, 37 

Mazzuca, L. M., Knapen, J. H, Veilleux, S., & Regan, M. W. 2008, 174, 337 
Mengel, S., Lehnert, M. D., Thatte, N., & Genzel, R. 2005, A&A, 443, 41 
Modigliani, A., Hummel, W., Abuter, R., et al. 2007, ArXiv Astrophysics e- 
prints 

Moorwood, A. F. M. & Oliva, E. 1990, A&A, 239, 78 

Mould, J. R., Huchra, J. P., Freedman, W. L., et al. 2000, ApJ, 529, 786 

Pantin, E., Lagage, P.-O., Claret, A., et al. 2005, The Messenger, 1 19, 25 

Peeters, E., Spoon, H. W. W., & Tielens, A. G. G. M. 2004, ApJ, 613, 986 

Piner, B. G, Stone, J. M., & Teuben, P. J. 1995, ApJ, 449, 508 

Quijada, M. A., Marx, C. T, Arendt, R. G, & Moseley, S. H. 2004, in 

Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) 

Conference, Vol. 5487, Optical, Infrared, and Millimeter Space Telescopes. 

Edited by Mather, John C. Proceedings of the SPIE, Volume 5487, pp. 244- 

252 (2004)., ed. J. C. Mather, 244-252 
Reach, W. T, Megeath, S. T, Cohen, M., et al. 2005, PASP, 1 17, 978 
Regan, M. W. & Teuben, P. 2003, ApJ, 582, 723 

Rosenberg, M. J. F, van der Werf, P. P., & Israel, F. P. 2012, A&A, 540, Al 16 
Sandage, A. & Tammann, G. A. 1981, Carnegie Inst, of Washington, Publ. 635, 
0,0 

Sanders, D. B., Mazzarella, J. M., Kim, D.-C, Surace, J. A., & Soifer, B. T. 

2003, AJ, 126, 1607 

Sarzi, M., Allard, E. L., Knapen, J. H, & Mazzuca, L. M. 2007, MNRAS, 380, 
949 

Schinnerer, E., Eckart, A., Quirrenbach, A., et al. 1997, ApJ, 488, 174 

Schwarz, M. P. 1981, ApJ, 247, 77 

Sersic, J. L. & Pastoriza, M. 1965, PASP, 77, 287 

Shlosman, I., Begelman, M. C, & Frank, J. 1990, Nature, 345, 679 



Siebenmorgen, R., Kriigel, E., & Spoon, H. W. W. 2004, A&A, 414, 123 
Smith, J. D. T., Draine, B. T., Dale, D. A., et al. 2007, ApJ, 656, 770 
Snijders, L., Kewley, L. J., & van der Werf, P. P. 2007, ApJ, 669, 269 
Telesco, C. M., Dressel, L. L., & Wolstencroft, R. D. 1993, ApJ, 414, 120 
Thornley, M. D., Schreiber, N. M. R, Lutz, D., et al. 2000, ApJ, 539, 641 
Tielens, A. G. G. M. & Hollenbach, D. 1985, ApJ, 291, 722 
Tully, R. B. & Fisher, J. R. 1988, Catalog of Nearby Galaxies (Catalog of 
Nearby Galaxies, by R. Brent Tully and J. Richard Fisher, pp. 224. ISBN 
0521352991. Cambridge, UK: Cambridge University Press, April 1988.) 
Uchida, K. I., Sellgren, K., & Werner, M. 1998, ApJ, 493, L109+ 
Vacca, W. D., Johnson, K. E., & Conti, P. S. 2002, AJ, 123, 772 
Vaisanen, P., Rajpaul, V, Zijlstra, A. A., Reunanen, J., & Kotilainen, J. 2011, 
ArXiv e-prints 

van der Werf, P. P., Snijders, L., Vermaas, L., Reunanen, J., & Hamelink, M. 

2006, The Messenger, 126, 11 
Verma, A., Lutz, D., Sturm, E., et al. 2003, A&A, 403, 829 
Whitmore, B. C, Chandar, R., Schweizer, F, et al. 2010, AJ, 140, 75 
Whitmore, B. C. & Zhang, Q. 2002, AJ, 124, 1418 



16 



